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Abstract 

In this paper we extend oblate and prolate Jaffe models into more general flat- 
tened Jaffe models. Since dynamical properties of oblate and prolate Jaffe Models 
have been studied by Jiang & Moss, they are not repeated here. 
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1 Introduction 

There is a long history of studying the structure of ellipsoidal galaxies by construction of 
self-consistent density-potential pairs. Some early profiles of modelling ellipsoidal galax- 
ies are spherical and purely empirical to fit the surface brightness of galaxies observed, 
for example, the de Vaucouleurs (119481) and Hubble (119301) profiles. The de Vaucouleurs 
profile can fit the brightness profiles of many giant ellipticals such as the giant El galaxy 
NGC 3379. In contrast, the data for NGC 4472 (a giant elliptical galaxy) cannot be fitted 
by the de Vaucouleurs profile. The Hubble profile fits the observed brightness profiles of 
some galaxies as well as the de Vaucouleurs laws. However, the Hubble profile predicts 
that the total mass of a galaxy is infinite. This must involve making some reasonable 
assumption, such as that the mass/light ratio is constant with position. Also, based on 
gravitational potential theory or Jeans's theorems, other spherical models have been con- 
structed by using potential-density pairs or described by using distribution functions, for 
example, the Plummer (119111) and King profiles. The Plummer profiles fit observations of 
globular clusters. Although the Plummer profiles can fit the brightness of some galaxies 
at finite R, they cannot fit the brightness profiles of galaxies at positions with large R, 
since the density distributions of galaxies typically decrease according to O(^) in regions 
of large r, where 2 < a < 4 (Binney Tremaine 1987 Mihalas fc Binney 1981). Ap- 



plying Jeans's theorems to spherical stellar systems, it follows that the potential-density 
pairs of Plummer's profiles can be reproduced by constructing a very simple distribu- 
tion function. The King profile is constructed from a distribution function (IMichie 19631: 
IMichie &: Bodenheimer 19631 King 1966 , I198ip and fits star counts in globular clusters 



and dwarf ellipticals (|Mihalas fc Binney 1981). The brightness profiles of some giant el- 



lipticals such as NGC 4472 can be also well fitted by the King profiles but none of the 
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King profiles can fit the observations for the giant El galaxy NGC 3379 as well as does 
the de Vaucouleurs law (Mihalas &: Binney 1981D, since the King profiles predict that the 



luminosity of a galaxy vanishes at some finite radius Rt- It is thus clear that neither 
empirical nor theoretical laws can describe all elliptical galaxies. 

Of course, there are also many other spherical models for the surface brightness pro- 
files of galaxies. Although spherical models can mimic the surface brightness of some 
galaxies observed, galaxies are not all spherical. Thus many axisymmetric non-spherical 
models have been constructed for our overall understanding of galaxy formation and evo- 
lution. The basic idea is to extend some spherical models to more general axisymmetric 
cases. The methods of this extension can roughly be classified into three groups. The 
first is to use the axisymmetric radius of \/ R 2 + (a + \z\) 2 in place of the spherical radius 
r = \AR 2 + z 2 , where a > 0. This was originally introduced by Kuzmin (119561) . It is 
necessary to discuss this device here although it is not directly relevant to ellipticals and 
it is used to model discs. The physical significance of this device is that at the point 
(R, — \z\) below Kuzmin's disc, the potential of Kuzmin's (119561) models is identical with 
that of a point mass located at distance a above the disc's centre. The second is to use the 
axisymmetric radius of \J R 2 + (z/q) 2 in place of the spherical radius, where q is the axial 



ratio, as in Binney's logarithmic models (Binney fc Tremaine 1987). The third is to use 



the axisymmetric radius of y R 2 + (V z 2 + c 2 + d) 2 in place of the spherical radius,where 
c and d are positive constants. This device was first introduced by Miyamoto and Nagai 
(I1975P for the Plummer potential. There are also other methods of constructing axisym- 
metric potential-density pairs. For example, a new potential-density pair can be obtained 
if a known potential is differentiated. Satoh's (I1980|) models is obtained by differentiating 
the potential of Plummer-Kuzimin's model n times with respect to c 2 if it is assumed that 



the axisymmetric radius of Plummer-Kuzimin's model is y R? + (V z 2 + c 2 + d) 2 as given 
above. 



Jaffe's (1983) spherical models can be flattened into the oblate models Piang 2000] ) and 



they can be also elongated into the prolate models Piang fc Moss 2002[) using Miyamoto 



and Nagai's device mentioned above. It have been known that the density of Jaffe's (119831) 
spherical model decays radially like r~ 4 at large distances and that, at large distances, the 
density of the oblate model given by Jiang (2000) decays radially like r~ 4 , except on the 
i?-axis, and like r -3 on the i?-axis. In contrast with the oblate model, the density of the 
prolate model piang Moss 2002) at large distances decays radially like r~ 4 , except on 



the z-axis, and like r~ 3 on the z-axis. The question addressed here is, can Jaffe's spherical 
model be extended into a model whose density at large distances decays radially like r -4 , 
except on the two principal axes, and like r~ 3 on the two axes? Yes. In this paper, we 
construct a class of more general flattened Jaffe models using a similar Miyamoto and 
Nagai's device, that is, replacing the spherical radius of the potentials of Jaffe's models 



with a more general axisymmetric radius of y (V-R 2 + o? + 6) 2 + (V ' z 2 + c 2 + d) 2 , where 
a, b, c and d are positive constants. Then we study the properties of the potential-density 
pairs of the general flattened Jaffe models in order to advance our overall understanding 
of galaxy formation and evolution. 
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2 Potential-density Pairs 

In order to answer the question mentioned in the last section, we recall Jaffe's spherical 
models with potential-density pairs as follows: 

GM r 

*(r) = —H-^—), (2-1) 
rj r + rj 

M r 4 

*> = ^mdv]?' (2 - 2) 

where and everywhere below, r = \J x 2 + y 2 + z 2 is one of three spherical coordinates 
(r,9,4>) which can be expressed by three Cartesian coordinates (x,y,z), M and rj are 
positive constants and G is the gravitational constant. Then we replace r in (12.11) by 



(VR 2 + a 2 + b) 2 + (Vz 2 + c 2 + d) 2 , 

where and everywhere below, a, b, c and d are positive constants, thus giving a family of 
general axisymmetric elliptical models with potentials 



GM / V VTO + b) 2 + (Vz^Tc 2 + d) 2 
$(i? 2 , z) = In ( / = | • (2.3) 

'(VR 2 + a 2 + b) 2 + (Vz 2 + c 2 + d) 2 + rj 



Of course, from ( 12. 3D . the densities 

M 

A^j r 4 (r + rj) 2 X 3 Y 3 



M rjAr* + r 2 Br 2 + (3r + 2r J )r J C 

p\tX ,Z) - \2V3V3 y^^) 



can be found, where 

A = bX 2 Y 3 + a 2 bY 3 + c 2 dX 3 , 
B = X 3 Y 3 + bX 2 Y 3 + a 2 bY 3 + c 2 dX 3 , 
C = a 2 X(X + b) 2 Y 3 + c 2 X 3 Y(Y + d) 2 , 
t = y/(X + b) 2 + (Y + d) 2 , X = VR 2 + a 2 , Y = Vz^+ 



2 



Obviously, p(R 2 ,z) > for all positive constants a,b,c,d and rj. It is also found from 
(12. 3p and (12.41) that the model (I2.3p - (I2.4I) degenerates into the oblate model as a and b 
limit to zero and the prolate model as c and d go to zero, respectively. It will be shown 
in the next section that the density of the model (I2.3p - (12.4p at large distances decays 
radially like r -4 , except on the two principal axes, and like r~ 3 on the two axes. 



4 



Z. Jiang, D. Fang, H. Liu, D. Moss 



3 Properties of Densities 

For the model (I2.3l) - (l2.4p . the density p(R 2 , z) at large distances has a radial dependence 
~ r~ 4 , except on the two principle axes, and like r~ 3 on the two axes, i.e., 

The ratio, a, of i?-axis to z-axis extent of the contours of finite p(R 2 , z) near the origin is 
of the form 

2 _ a 2 [3ad,TQ (r + rj) + c(c + d)(3ad + 26c + 5ac)(3r + 2rj)] , . 

c 2 [46cr 2 (r + o) + a(a + 6) (ad + 46c + 5ac)(3r + 2tj)\ ' 

where ro = y(a~+ 6) 2 + (c + d) 2 . It is worth mentioning that the ratios in the oblate and 
prolate Jaffe models are degenerate forms of (I3.3p . In fact, by substituting 6 = into 
( 13. 31) and letting a go to zero, ( 13. 3D becomes the ratio in the oblate Jaffe models; similarly, 
( 13. 3j) with d = can be changed into the ratio in the prolate Jaffe models as c tends to 0. 
It is also below found by element integration that the total mass of the general flattened 
Jaffe models is M. 

In order to consider the total mass of the more general models (I2.3l) - (l2.4p with positive 
constants a, 6, c and d briefly, it is necessary to introduce the following notations: 

r _ ^ rjcH rj(3T + 2rj) c\Y + df 

1 ~ r(r + rj) 2 r(r + rj)Y 3 r 3 (r + o) 2 1 Y 2 h { ' 

r,jb rja 2 b 

h 



(r + rj) 2 (X + 6) T(T + rj)X 2 (X + 6) 



r.M 2 d + 



t(t + rj)(X + b)Y 3 



1 

+ 



a 2 6 bc 2 (Y + df 



X (X + b)Y 2 



(3.5) 



r 3 rjT 2 rj(r + rj) 2 
where X, Y and r are the same as in (12.41) . thus, by (12.41) . giving an identity 

Anr jP {R 2 ,z)XT 

M(X + b) =h + h - (3 ' 6) 

Integrating (13.41) about r from \/{a + b) 2 + (Y + d) 2 to +oo, it follows that 

Y 3 + c 2 d , / v/(a + 6) 2 + (Y + d) 2 + rj 
i\dr = — — In — ====== 

V^FTOW ^ 3 V y/(a + b) 2 + (Y + d) 2 

rj[(Y + d) 2 (F 2 - c 2 ) + (6 2 + 2a6)F 2 ] 
F 2 [(a + 6) 2 + (Y + d) 2 ] [ v/(a + 6) 2 + (FTd) 2 + r j] ' 



(3.7) 
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Integrating (13.51) about r gives 



hdr 



c 2 (Y + d) 



21 



r 2 Y 2 



rjb(X + b) 



((Y + d) 2 -rj)(r + rj) 



1 - 



c 2 (Y + df 



rjY 2 



rjb 



bc 2 d 



(Y + d) 2 - rj Y 3 
rja 2 6 6c 2 (X + 6) 6c 2 (X + 6) 



mr, Y + d) — — arccos 



r 2 (r + rj)X r 2 Y 2 
for Y + d ^ rj, where and everywhere below, 



+ 



rjrY< 



constant 



g{r,Y + d) = < 



^(Y+d) 2 -r 2 



Jr 2 ,-(Y+d) 2 1 I rj+Y+d [7 
V J \V rj-Y-d->jT 

With the help of (E3D dSIHD, 



™Y + d> rj , 

( j rj+Y+d / T -Y-d \ 
I V rj-Y-d^y T + Y + d I 

: In I -^== -= I as Y + d < r j. 



+Y+d 



p{R 2 ,z)RdR 



A7crjp{R 2 ,z)XT 



y/(a+b) 2 +{Y+d) 2 M(X + b) 



dr 



Y 3 + c 2 d , / v/(a + 6) 2 + (Y + d) 2 + r j 
In 1 



Y 3 



V(a + b) 2 + (Y + dp 



rj[(Y + d) 2 (Y 2 - c 2 ) + (6 2 + ab)Y 2 ] 



Y 2 [(a + 6) 2 + (Y + d) 2 ] + &) 2 + (Y + d) 2 + o] 
c 2 (Y + d) 2 ] rj6 



1 - 



+ 

c 2 (Y + d) 2 
rjY 2 

c 2 (Y + d) 2 ' 
rjY 2 



1 - 



rjY 2 J (Y + d) 2 -r 
bc 2 d 



rjb 



(Y + d) 2 - rj Y 3 



I {Y + d,rj) + 



be 2 



rjY 2 



rjb(a + b) 



+ 



c 2 (Y + df 



((y + d)2 _ r 2 )(v/(a + 6) 2 + (Y + d) 2 + rj) 
6c 2 d 



r 2 6 



rjY 2 J (Y + d) 2 - r 2 Y 3 



(3i 



(3.9) 



be 2 (a + b) 



be 2 

-arcsm 



g{^{a + b) 2 + (Y + d) 2 , Y + d) 
Y + d 



r J Y 2 v ^T6) 2 T(YTd) 2 ^ 3 

6c 2 (a + b) 



v/(a + 6) 2 + (Y + d) 2 



+ 



((a + 6) 2 + (Y + d) 2 )Y 2 



(3.10) 
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implies that the total mass is finite for the model given by f)2.4p . as follows. It can be 
deduced from (I3.10p that 



[ ( r z) mr) dz = + m ( ^q^+O^+r; 



fcs(Y + d) 



y 



fez 

— arcsm 



/ (y + d, rj) - g(y/(a + b)* + (Y + d) 2 , Y + d) 
Y + d 



constant 



^(a + b) 2 + (Y + d)\ 
for any positive c, since the following equalities hold: 

| Li i ^-(^™ + ^ 

^ \ y \ V(y + d) 2 + (a + 6) 2 
y 3 + c 2 d ^ / ^/(a + b) 2 + (Y + d) 2 + r / 

" ^ 3 11 v V( a + & ) 2 + + rf ) 2 

r,j(Y + d) 2 (Y 2 - c 2 ) 
~ Y 2 [(a + b) 2 + (y + d) 2 ] &) 2 + (y + dj 1 + O] ' 



and 
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d 



— 6 2 (i + -)/o(y + rf,o) 



y 



6c 2 (Y + d) 3 - r 2 fe(Y 3 + c 2 d) r 

y 3 [ (F + ci )2_ r ,2 ] " 



(Y + d,rj) + 



rjbiY 2 - c 2 " 



Y 2 [(Y + d) 2 - r 2 ]' 



d 



fe(i + -)^( v / (a + 6) 2 + (y + d) 2 , y + d) 



6c 2 (y + d) 3 - rjb(Y 3 + c 2 d) 



(^(a + b) 2 + (Y + d) 2 , y + d) 



y 3 [(y + d) 2 - r 2 ] 

rj6(a + 6)(y 2 - c 2 ) 
+ y 2 [v^+6) 2 + (y + d) 2 + rj}[{Y + d) 2 - r 2 ] 

fe(a + fe)(Y 2 -c 2 ) 

yy(a + 6) 2 + (y + d) V(a + fe) 2 + (y + ^) 2 + o] 



d \bz 
dz I y 



arcsm 



y + d 



6c 2 



y + d 



arcsm 



V / (a + 6) 2 + (y + d) 2 ; 

&(a + 6)(Y 2 -c 2 ) 



+ 



V(a + 6) 2 + (y + d) 2 / y 2 [(a + &) 2 + (y + d) 2 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



(3.15) 
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Furthermore, it can also be shown from (13.111) that 

poo roo 

4tt / / p(R 2 , z)RdRdz = M, 
Jo Jo 

i.e., the total mass is still M for the more general models fl2.3l) - fl2.4|) . 

Remark: The potential Q(R 2 ,z) ~ —GM/r + 0(l/r 2 ) at large distances and the 
density is obtained from the potential by Poisson's equation AQ(R 2 ,z) = A7cGp(R 2 , z). 
The divergence theorem then guarantees that the volume integral of the analytic p(R 2 , z) 
gives a total mass of M, i.e., the Cauchy value of the volume integral can be easily obtained 
by using the divergence theorem and it limits to M as r — > +00. This idea is very simple. 
In spite of this, it is necessary to do a set of integrations here since the method used to do 
these integrations is similar to that of deriving other formulae of the edge-on projected 
surface densities for the flattened Jaffe models and it is very useful to the readers for their 
understanding of how to obtain those results. 
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